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Hamiltonian lattice gauge theory: wavefunctions on large lattices 

J. B. Bronzan'** 

'^Department of Physics and Astronomy, Rutgers University 
P.O. Box 849, Piscataway, New Jersey 08855-0849, U.S.A. 

We discuss an algorithm for the approximate solution of Schrodinger's equation for lattice gauge theory, using 
lattice SU(3) as an example. A basis is generated by repeatedly applying an effective Hamiltonian to a "starting 
state." The resulting basis has a cluster decomposition and long-range correlations. One such basis has about 
10* states on a 10 X 10 X 10 lattice. The Hamiltonian matrix on the basis is sparse, and the elements can be 
calculated rapidly. The lowest eigenstates of the system are readily calculable. 



The approximate solution of Schrodinger's 
equation for lattice gauge theory presents a num- 
ber of problems, especially when it is to be carried 
out on a large lattice. We discuss these issues for 
pure SU(3) field theory, sometimes using the ex- 
ample ofalOxlOxlO lattice, where there are 
3000 link degrees of freedom. The addition of 
matter degrees of freedom is straightforward. 

1. CHOICE OF THE HAMILTONIAN 

The computations of SU(3) wavefunctions is 
based on the Kogut-Susskind Hamiltonian for 
SU(3) lattice gauge theory. |l[] 
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This Hamiltonian commutes with local gauge 
transfromations, and it is only its locally gauge- 
invariant eigenstates that are physically signifi- 
cant. One way to eliminate states having color 
charge on the lattice is to fix the gauge. The 
resulting Hamiltonian has fewer degrees of free- 
dom, but is very messy. ||^ For example, the elec- 
tric term in the Hamiltonian becomes 



= ^ {PlJL,a{yal.bl,cl)JL,a{ya2,b2,c2) 

+ PlJL.aiZal, hi, cl)JL,a{Za2,b2,c2) 



+(33'JR,a{ZalM,cl)>JB.,a{Za2,b2,c2) 
+P4 [<JL,a{yalM,cl) + >JL,a{Zal,bl.cl) 
X [JL.a{ya2,b2,c2) + JL,a{Za2.b2.c2)] 



+(^5 [<jR,a{yal,bl,cl) 
X [JB.,a{ya2.b2,c2) + 
+/36 [<JF/.,a{Zal,bl,cl) 
X>7l, 0(202, 62, c2)}- 



+ <jR,a{Zal,bl-cl)] 
<jR,a{Za2,b2,c2)] 
— •JL,a{ZalM.cl)] 



(2) 
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The functions /3 are nonzero for degrees of free- 
dom located arbitrarily far from one another. 
Here is one of them: 

/34(ai, 61, ci; a2, 62, C2) = min{cl, c2). (3) 

The locality and homogeneity of the Kogut 
Susskind Hamiltonian Eq. |^ have been lost in 
Eq. |2[ However, these characteristics are cru- 
cial to our construction of basis states, and we 
therefore opt to use Eq. |^. This implies that the 
states we construct will be contaminated with 
color. Our strategy is to construct eigenstates 
of Eq. |l| for a range of expectation values of the 
lattice color operator = J2iQi- (The sum 
is over sites i, where Qi^a a-i'e the generators of 
gauge transformations.) We then extrapolate our 
results to color zero. 

2. BASIS STATES 

Basis states are products of factors or 0q, 
for each degree of freedom. These factors are 
"harmonic oscillator" states on the SU(3) man- 
ifold.]^ They transform as scalars and octets, re- 
spectively, under SU(3) transformations. Higher 
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tensor states exist, but have proven unimportant 
in work on small lattices. We omit them at this 
stage of our work. 

The construction of a basis begins with a 
"starting state" |0) in which each link has the 
same factor (j). This state has a parameter in, the 
width of the harmonic oscillator Gaussian, and it 
is chosen to minimize the expectation of H at the 
couphng we choose. At strong couphng the state 
is the true ground state. At weak coupling the 
expectation value is 26% above the ground state 
energy. At weak and intermediate couplings the 
starting state has nonzero color density, which 
must be eliminated later by extrapolation. 

It is not feasible to add basis states by using 
the nine states t^, (f>a as a basis on each link. On 
alOx 10x10 lattice there are 3000 links and there- 
fore 9^°°" lattice states are generated in this man- 
ner. It is better to set up a basis using dynam- 
ical information by making the following query: 
If I?') is a lattice state, what is the optimal state 
A\^) to add so that a linear combination can be 
formed having given color (Q^) and minimum en- 
ergy? The (approximate) variational choice for A 



A = H + XQ^ 



(4) 



Here A is a parameter. As we tune it, we change 
both (Q^) and the variational energy. The rest of 
the basis is generated by repeated application of 
A to the starting state: 



= (^)"|0). 

3. EXTRAPOLATION 
COLOR 
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Using basis |^, it is possible to make a simple 
test of the extrapolation procedure to zero color. 
We do this at weak coupling where perturbation 
theory provides a comparison benchmark. The 
basis we use is the smallest one that allows study 
of the problem: $o and <I>i. It is convenient to 
write the color and energy in terms of variables x 
and y. 



4L3 16^3/2 
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L is the number of links on the edge of the lat- 
tice. Variables a;, y fall on a tilted eccentric ellipse 
{h/a = 0.41), given by the equation 

8a;2 - IQxy -f 2\y^ = 13. (7) 

As A sweeps from — oo to +oo, x varies from 
+ y/2l/^ to -v/21/8, and y moves over the lower 
branch of the ellipse from + ■\/8/21 to — ^8/21. 
The simplest extrapolation is to use a straight 
line passing through these endpoints to reach 
(Q^) — 0. When this is done, the extrapolated 
ground state energy extrapolated to the physical 
sector becomes Eq = which is 8% above 

the exact weak-coupling result i?o = 19. IL^. This 
is to be compared with the 26% excess when only 
starting state $0 is present. (This state has color 
density {Q^) / ^ 48 / .) 

4. THE CLUSTER EXPANSION 

The Kogut-Susskind Hamiltonian is a sum of 
operators that are either local or nearly local, and 
the same is true of the color operator Q^. We 
denote these terms generically as Ai; then Eq. ^ 
becomes 



|$„) = A^iA^2---A^„\0). 



(8) 



Terms involving different links commute, so we 
can rearrange this into the form 



|<I>„) = [A.aAa ■ ■ ■][A,c,A 



|0), (9) 



where the indices are a permutation of il, . . . , in, 
and the factors within a given bracket share a 
degree of freedom with another factor within the 
same bracket, but not with any factor in other 
brackets. These bracketed factors produce clus- 
ters of excitation of the harmonic oscillator state 
^ 4>a on the lattice. The cluster representation 
is a useful framework for constructing a basis be- 
cause A is a sum of nearly local and translation 
invariant operators. Clusters are therefore small 
(locality) and identical wherever they appear on 
the lattice (translation invariance) . Our approach 
is to compute and store cluster amplitudes, and 
then use an elaboration of Eq. ^ to compute lat- 
tice wavefunctions from them. This approach is 
not available with the gauge fixed Hamiltonian |^. 
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There is one further deomposition of <&„ that 
we choose to make. In Eq. ^ different index as- 
signments lead to different numbers of clusters. 
We separately group the terms having Uc clusters 
into a substate ^n,n^ and retain these different 
linear combinations as separate cluster states. 

When a single term in operator A is applied 
to the starting state, it can have several effects, 
including the creation of a state with a pair of 
excitations having contracted color indices. 



A,|0) =^a(/)(l)---(/)„(j)---0„(A:) 



(10) 



where j and k designate links appearing in Ai. 
Application of a second term in A can either 
change the first pair, or create a second cluster. 
A particularly important change occurs when the 
second application extends the distance between 
links j and k to build longer-range correlations. 
This is a contribution to $2,1 and we know that 
such terms are crucial in the scaling region. Un- 
fortunately, such long-range terms are apt to be 
swamped in the matrix elements. For example, in 
the norm of the state $2, the contribution of $2,2 
dominates the contribution of $2,1 by a factor of 
the lattice volume because the two clusters can 
range independently over the lattice. 

We therefore make adjustments to the basis im- 
plied by Eq. ^: 

1) As mentioned, we separate $„ into distinct 
states $ri,nc to evade swamping. 

2) In each cluster, we retain only a one-pair ex- 
citation like that discussed. When this is done, 
clusters are the same thing as pairs. Pairs are 
computed iteratively. That is, to obtain the pairs 
of order N, we apply A to the clusters of or- 
der A'^ — 1, and store amplitudes. The number 
of distinct kinds of pairs of clusters can be no 
more than 3L{L + 1)^/2; this number of clus- 
ters presents no computational problem on a 
10 X 10 X 10 lattice. 

The basis we construct in this manner is a 
stripped-down version of Eq. H, but it still is 
constructed dynamically, and contains long-range 
correlations. We believe these features are the es- 
sential ones for correct scaling. In addition, we 
believe it is important that our starting state is 
reasonably accurate at weak and strong coupling. 



We emphasize that the lattice states we con- 
struct have many pair excitations. For example, 
the state $137,111 has 111 pairs on the lattice. A 
related point is that when ric = 3L{L+ 1)^/2, ev- 
ery link on the lattice is excited. We include such 
states in our basis because the only way we can 
shift the ground state energy density is to tile the 
lattice with excitations. 

We anticipate including cluster states having 



n < 3L{L + iy/2 and n - r 
basis therefore grows like L^. 



^ N < L. The 



5. COMPUTABILITY OF MATRIX EL- 
EMENTS 

The lattice state ^n,na is a sum over states in 
which He pairs are placed on the lattice in var- 
ious ways. It is not immediately apparent that 
it is feasible to compute matrix elements of these 
states because of the enormous number of terms 
to be considered. Two features are of help here, 
and they can be discussed using the case of the 
computation of state norms. 

1) We saw in Section 4 that phase space favors 
those terms in which a maximum number of pairs 
can be positioned independently on the lattice. 
In the case of the norm, contributions in which 
the set of pairs in the bra is identical to the set 
of pairs in the ket dominate other contributions 
by at least one power of the lattice volume. We 
retain only these leading contributions. 

2) There is still danger that the norm will be 
uncomputable. It has the stucture 



pairs 



(11) 



Here i enumerates pair type, Ni is the number 
of applications of A required to produce pair i, 
and there are Uc pair types in each contribution. 
a({i}) is an amplitude assembled from the pair 
amplitudes Xi computed by the iterative proce- 
dure outlined in Section 4. The number of pair 
types is of order i^, so the number of terms in 
Eq. is the same as the number of partitions of 
n into non-negative integers, or 



(L3 ~ l)!n!' 



(12) 



4 



This is an impossibly large number on a lattice of 
any size. Fortunately, the amplitudes a in Eq. |ll| 
depend on the pair amplitudes in a manner that 
allows us to use the distributive laws of arithmetic 
to sum before multiplying the pair amplitudes. 
This possibility is illustrated by the choice 

a^X{{x.r^/{N.)\ (13) 

i 

Then 




Similar simplifications occur in the computation 
of other matrix elements of A. These simplifica- 
tions are discovered by writing generating func- 
tions for matrix elements. Manipulation of these 
matrix elements teach us how to interchange mul- 
tiplication and summation of pair amplitudes. 
Preliminary numerical work onalOxlOxlO lat- 
tice indicates that the computation of the Hamil- 
tonian matrix is not taxing. 

6. CONCLUSIONS 

We are working at present onalOxlOxlO 
lattice, in which n^. ranges up to 1,500, so that 
the last states tile the entire lattice with pairs. 
We allow n to exceed ric by as much as 6, so that 
pairs separated by up to 6 lattice spacings appear 
on the lattice. There are thus 9001 states in our 
basis. The Hamiltonian has ~ 4 x 10^ distinct 
matrix elements, but of these only 162,000 are 
nonzero. On a workstation nonzero matrix ele- 
ments are computed in less than one minute. We 
have previously discussed the problem of finding 
the low-lying states of such a matrix[Q. 
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